Fundamental Limits of Electronic Nonlinear-Optical Phenomena 
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Using sum rules and a new dipole-free sum-over-states expression, we introduce a method for 
calculating the fundamental limits of the dispersion of the real and imaginary parts of all electronic 
nonlinear-optical susceptibilities. As such, these general results can be used to study any nonlinear 
optical phenomena at any wavelength, making it possible to push both applications and our under- 
standing of such processes to the limits. These results reveal the ultimate constraints imposed by 
nature on our ability to control and use light. 

PACS numbers: 42.65.An, 33.15.Kr, 11.55.Hx, 32.70.Cs 
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I. INTRODUCTION 

The interaction of light with matter is of funda- 
mental importance in studying materials; and, lies 
at the heart of many critical technologies that span 
telecommunications, optical data storage, three- 
dimensional nano-photolithography,0 0| and making 
new materialsloj for novel cancer therapies. Be- 
cause the strength of interaction, as quantified by the 
nonlinear-optical susceptibility, governs the efficiency of 
an application - and indeed whether or not a particu- 
lar application is practical, making materials with ever 
larger nonlinear susceptibility has been the central focus 
of research in organic nonlinear optics. Is there a funda- 
mental limit to the susceptibility? Our work shows that 
nature imposes such an upper bound on all optical phe- 
nomena. Prior work determined the fundamental limit of 
only the off-resonant susceptibility. H, 0,0,0,0,01 In 
the present studies, we calculate the most general case: 
the fundamental limits of the dispersion of the real and 
imaginary parts of all electronic nonlinear-optical suscep- 
tibilities, which can be used to study or apply any non- 
linear optical phenomena at any wavelength. Our work 
provides an understanding of the constraints imposed by 
nature on our ability to control and use light, and lays 
the foundation for developing better materials and novel 
applications. 

We focus on the second-order nonlinear-optical suscep- 
tibility of a molecule, often called the hyperpolarizability 
(3, as an example of the process we use for finding the 
fundamental limit of any nonlinear susceptibility. Fur- 
thermore, since our goal is to set an upper bound, we will 
only consider the largest tensor component, Pxxx- The 
electronic nonlinear-optical susceptibilities are calculated 
using perturbation theory in the dipole approximation, 
which yields an expression that incudes the excited state 
properties of all the quantum states of the system. This 
theoretical result is called a sum-over-states (SOS) ex- 
pression, and for (3 is given bv:[l3| 
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where — e is the electron charge, Xnm the n, m matrix el- 
ement of the position operator, Axno = Xnn — xqq is the 
difference in the expectation value of the electron posi- 
tion between state n and the ground state, D~^{uJa,ujf3) 
gives the dispersion of P (defined later) and htUi, are the 
photon frequencies. The primes indicate that the ground 
state is excluded from the sum and the permutation op- 
erator Puj^.uji^ directs us to sum over all six frequency 
permutations. Since the dipole moment of the molecule 
is proportional to the position (px — —ex), we loosely 
call Xnrn the transition moment and Xnn the excited state 
dipole moment. The first and second terms in Equation 
Hare called the dipole and the octupolar terms. 

Equation ^ is a function of an infinite number of ma- 
terial parameters, Xnm and Eno, so the maximum value 
of (3 for each possible pair of photon energies requires an 
optimal set of transition moments and energies. The sum 
rules, which are directly derivable from the Schrodinger 
Equation, are relationships between the transition mo- 
ments and energies. In the present work, we apply the 
sum rules to put the SOS expression of the nonlinear op- 
tical susceptibilities into a form that can be maximized 
to calculate the ultimate nonlinear-optical susceptibility. 



II. THEORY 

The generalized Thomas-Kuhn sum rules, derived di- 
rectly from the Schrodinger Equation without any ap- 
proximations, relate the matrix elements and energies to 
each other according to:@ 
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where m is the mass of the electron, and N the number 
of electrons. The sum, indexed by n, is over all states of 
the system. Equation |21 represents an infinite number of 
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equations, one for each value of m and p. As such, we 
refer to a particular equation using the notation (m^p). 

Defining Eij = Ei — Ej , we eliminate the dipole term 
using Equation {m,p) with m 7^ p: 
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Equation 13 can be rewritten by explicitly expressing the 
n = m and n = p terms, setting p = 0, and multiplying 
both sides by a:om:Ii3| 
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Substituting Equation^ into Equation we get. 
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+ LUi LU2 for the three previous terms] , 

and where oj^g — ^mg ~ *r„/2. huj^g is the energy dif- 
ference between states Em and the ground state and 
is the damping width. Equation|Slis called the dipole-free 
expression or the reduced hyperpolarizability, 15j where 
the second term in brackets implicitly accounts for the 
dispersion of all dipolar terms. 

We use the ansatz that at most two states contribute to 
a molecule's nonlinear susceptibility near its fundamental 
limit. [H, [13 (We str ess that a three-level SOS model 
does not demand that the sum rules are truncated if the 
energies of the higher-lying state are high enough, so no 
truncation pathologies result. Details will be described 
in a future publication.) Using Equation [31 
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where D^^{uji,uj2) is the term in brackets. 

To further simplify Equation[71so that we can calculate 
the extrema, we again use the sum rules to express the 
matrix element product /ioi/ii2/^20 hi terms of /ioi. This 
leads to the exact result (without truncating the sum 
rules) is given by 
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where E — Eiq / E20 and max is the maximum allowed 
transition dipole moment to the first excited state, which 
is related to the number of electrons, iV, and the energy 
to the first excited state, i?io, according to 
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Having exhausted all the sum rules that meet the 
three criteria for physical soundness. |l6l IT^ we assume 
that the remaining energies and transition moments are 
independent, so /3 (via Equation (HJ is maximal when 
fJ-io = \^IJ'MAX- Combining Equations |H1 and El 
with fiio = '^VSfiMAX, we get 

Pxxx{cOl,U:2) < P^'^"" ■ \^=^ ■ D''^{UJ^,U2) (10) 

where Pq^^-^ is the fundamental limit of the off-resonant 
^hyperpolarizability: |lflj 
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The maximum value of Pxxx for any second-order process 
at any wavelengths uji and 102 is given by Equation 1101 - 
a function of E, Eiq, r„o, and N. 



III. DISCUSSION 

All second-order nonlinear optical processes are 
bounded by Equation ^1 if the three-level ansatz is 
obeyed. Figure ^ shows the real and imaginary parts of 
the maximum allowed hyperpolarizability as a function 
of the energy of the two incident photons, normalized to 
the maximum off- resonant value Pq'^^'^ . Energies of the 
two dominant states used in these plots are Eiq = leV, 



E20 = 2eV and Vf^H^ ^ y^^^"^' = 100 meV (a 
100 meV width is a common approximation for organic 
molecules) . 

There are two dominant peaks in the fundamental limit 
of the real part of Pxxx- The electrooptic hyperpolariz- 
ability, (3xxx{0,u!), peaks at just under 80 times the off- 
resonant fundamental limit while the second harmonic 
generation hyperpolarizability, Pxxxi^,^), peaks at just 
under 50 times the off-resonant limit. Using these ex- 
cited state energies and damping factors, the ultimate 
hyperpolarizabilities for an A^-electron system is 
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FIG. 1: The dispersion of (a) the real part and (b) the imag- 
inary part of fundamental limit of /3 normalized to the off- 
resonant limit with Ew = leV , E20 = 2eV and rfo^^" = 
rfcT'^*^ = WOmeV for a one-electron system. 




FIG. 2: (a) The real part of fundamental limit of /3 normalized 
to the off-resonant limit with _Eio — leV as a function of 
£20, for various values of rfo^^*^ = FIo^^"; and (b) The 
real part of the fundamental limit of the second harmonic 
/3 normalized to the off-resonant limit with £10 = leV as a 
function of £20, for various off-resonant pump wavelengths. 



For the off-resonant fundamental limit of /3(0, 0), the 
second dominant excited state energy is required to be 
infinite. Since the energy level spacing typically gets 
smaller for higher energy states, most molecules should 
fall far short of the limit. Indeed, this could explain the 
factor of 30 gap between the best molecules and the off- 
resonant hmit.jTslIT^ 

Figurc|3i shows a plot of the real part of the fundamen- 
tal limit of the resonant electrooptic hyperpolarizability 
/3(0, huj = Eio) normalized to the off-resonant limit (with 
EiQ = leV) as a function of E20, for various values of 
YFWHM ^ pFjyffA/^ expected, the maximum reso- 
nant value increases as the width decreases. However, 
the hyperpolarizability peaks when the second dominant 
excited state is close in energy to the first dominant ex- 
cited state. In the limit of zero width, (3 peaks when the 
two dominant states are degenerate. Such excited state 
energy spacing is common in organic molecules, so it may 
be possible to identify molecules with an ultra-large res- 
onant response using this strategy. 



The structure required of a molecule to be at the fun- 
damental limit is thus different on and off resonance. 
Our present calculations, in the zero-frequency limit of 
are summarized by the series of plots in Figure 
Eb, which show the fundamental limits for several pho- 
ton energies between and OAeV with a damping factor 
of F„o = lOOmeV. For huj — 0, the fundamental limit 
clearly approaches the normalized value of unity, as pre- 
viously calculated. For each off-resonant photon energy, 
the fundamental limit increases as tku — > Eiq/2 and the 
second excited state energy approaches infinity. 

The limits in the dispersion of the real and imaginary 
parts of the third-order susceptibility, 7(0;!, a;2, a;2), can 
be calculated in a similar way using the three-level ansatz 
in dipole-free form. [Tsjl The sum rules can then be used 
to re-express the results in a form analogous to Equations 
[TUl A full discussion of third-order susceptibilities will be 
presented in a future paper. 

It is interesting to apply our results to typical 
molecules. For example, consider a molecule with 50 
conjugated electrons whose first excited state energy is 
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Eio = 1.55 (800 nm) and T = 100 meV. According 
to Equation Illl fT^ the ofF-resonant value of the hyper- 
polarizabihty yields (3^^^^ = 90, 500 x 10^^° esu. Using 
our new theoretical results, the resonant electrooptic hy- 
perpolarizability is I3max{Q,1T'UJ — Eiq) — 7,240,000 x 
lO^'^^esu. For a polar-aligned bulk material made of such 
molecules, at a concentration of 10^^ cm~'^, the funda- 
mental limit of the ofF-resonant electrooptic coefficient is 
rMAx{0,0) = 3,000pm/V and the resonant case yields 
rMAx{0, Tt-^ = Eio) = 240, OOOpm/y. Higher concentra- 
tion systems, or ones with narrower widths, would have 
even higher values. As such, our calculations show that 
there is room for considerable improvement over today's 
best materials. 

Our theory is general and few approximations have 
been used in calculating the dispersion of the limits. We 
start with a dipole-free form of the SOS expression for 
f3 and apply the ansatz that the oscillator strength is 
dominated by two excited states when the susceptibility 
approaches the fundamental limit. So in our calculation, 
the three-level model is not an approximation, but an ex- 
act model for a system at the limit. 16, 17] Furthermore, 
since this does not demand the sum rules to be truncated 
to three levels, no truncation pathologies result. The fact 
that most real molecules are not well approximated by 
a three-level model may explain the factor-of-thirty gap 
between the fundamental limit and the best molecules. 
The fact that a clipped harmonic oscillator yields a value 
of (3 that is near the fundamental limit 19] shows that the 
calculated limits are not large overestimates. The set of 
all theoretical calculations and measurements are consis- 
tently below the fundamental limit, leading credence to 
the theory. 

We only consider processes that are described by a 
sum-over-states expression which inherently assumes the 
electric dipole approximation (i.e. the wavelength of light 



is small compared with the size of the quantum system 
under study). While it is often the case that higher-order 
terms or magnetic moments may be larger than lower- 
order terms, the fundamental limits of the higher-order 
terms should not be larger than the fundamental limits of 
lower-order terms. Secondly, we only consider the largest 
diagonal components of /? and 7. Short of unexpected 
new Physics, we would not expect that the maximum 
limit of an off-diagonal component of a susceptibility to 
exceed a diagonal one. Note that our calculations include 
resonance, so Kleinman Symmetry is expected to be bro- 
ken, leading to potentially large off-diagonal components 
of susceptibilities that are disallowed off resonance. To 
treat such cases, one would need to use off-diagonal sum 
rules. 

In conclusion, we have applied the sum rules to deter- 
mine the fundamental limits of the dispersion of any com- 
plex susceptibility, thus being applicable to any electronic 
nonlinear-optical process. We find that to attain the fun- 
damental limit, off-resonant processes require a different 
molecular structure than on resonance, so one approach 
of molecular engineering does not suit all applications. 
Our work can be applied to calculating and optimizing 
the figure of merit of materials for specific applications; 
provides a method for comparing measured values of /3 
and 7 without having to resort to determining Pq and 70, 
which requires the use of unreliable dispersion models; 
and gives guidance for designing new materials and how 
their dispersion can be tuned for for a desired response. 
Clearly, there is still room for substantial improvements 
in all nonlinear molecular susceptibilities, especially on 
resonance in systems with closely-spaced excited states, 
such as provided by octupolar molecules. |20( 
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